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("■^ ' We comment on some open questions and theoretical peculiarities in Tsallis nonextensive sta- 

f^ , tistical mechanics. It is shown that the theoretical basis of the successful Tsallis' generalized ex- 

^S| . ponential distribution shows some worrying properties with the conventional normalization and the 

escort probability. These theoretical difficulties may be avoided by introducing an so called incom- 
plete normalization allowing to deduce Tsallis' generalized distribution in a more convincing and 






fT , ' consistent way. 
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I. INTRODUCTION 

It is well known that Boltzmann-Gibbs statistics (BGS) is inadequate for treating some complex 
systems. These are systems with complex or long term interactions and correlations, systems showing 
often distribution laws different from the usual ones (Gauss, Poisson), systems in chaotic or fractal 
states and often related to nonextensive phenomena in energy, entropy, heat, and other quantities. 
Some examples of the failures of BGS are given in reference [|l| . We see that we need a new statistical 
mechanics fundamentally different from BGS. We can even conjecture that a new kind of statistical 
theory may be necessary for complex random phenomena because the validity of the actual statistical 
method is subject to some ideal conditions [|2[. A brief discussion on this topic will be given in the 
present paper. 

In 1988, in a historical paper [^Q, Tsallis founded a nonextensive statistical mechanics which, in 
its most recent version, gives following canonical distribution functions : 



[l-(l-g)yV(e.-f/)]| 

P^ = H? 



1 

-9 



with 



Z = t[^~il-q)^ie,-U)f^ (2) 

i 2-ii Pi 

where [y]-^ = y if y > and [y]^ = otherwise (Tsallis cut-off), i is a state point in phase space, 
w the number of all accessible phase space points for the considered system, Cj the energy of the 
system at the state i, (3 the Lagrange multiplier of the constraint on the internal energy U given by 

which is used by Tsallis to overcome some theoretical difficulties @,|^. Eqs. (|l]) and (Q) mean 

w 

T.P^ = h (4) 

i 

which is logical if we consider the fact that w represents all possible states of the system. It should 
be noticed that, in the limit q = I, Eq. (|T]) becomes BGS, i.e. Pi = -|exp(— /?ej). 

We refer to the generalized exponential distribution Eq. (|l]) as Tsallis distribution function (TDF) 
which is proved to be indeed very useful and efficacious in describing some systems with complex 
or long term interactions. Many successful and convincing applications was published over the last 
10 years concerning different systems showing peculiar distribution laws (the reader is referred to 
references |]l|,||J3 for updated comments on this subject). 

In spite of the successes oiTDF, researchers continue to work on its theoretical foundation, because 
there are still various open questions. One notes peculiar theoretical properties which sometimes are 
not very easy to understand and deserve to be investigated. In this paper, we would like to discuss 
some fundamental aspects of the last version of Tsallis theory and present some observations. We also 
show a possible solution to the problems with a modifications in the theoretical basis of nonextensive 
statistics. 



II. BRIEF HISTORY OF TDF 

Let us begin by a brief review of the historical path of TDF. 

In 1988, Tsalhs proposed the following generalized statistical entropy : 

S = -k ^^'='^\ {qeR) (5) 

l-q 

which he maximized according to Jaynes principle with two constraints : Eq. @ and the mathe- 
matically suitable expectation value of energy U : 

w 

U = J2Pi(^i- (6) 

i 

This definition is completely legitimate if pi is considered as probability. This approach led to the 
first version of TDF, i.e. 



Vi = ^ — (7) 



[l-(g-l)/?e-i-^ 
Pi = 

with 
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Z = J2[l-iq-l)Pe,]f' (8) 

i 

Later, it was found that Tsallis — 1 had some shortcomings such as the thermodynamic instability 
of entropy and the incapacity to deduce some power-law distributions P,p|-p!l|]. But quite before the 
first published criticisms P,pT]] on this theory, another version of TDF was proposed by Curado and 



Tsallis |T^. This second version, say, Curado-Tsallis, replaced Eq.(^ by 

w 

u = Y.vW (9) 

i 

Curado-Tsallis formalism allows to get a mathematical elegance with the conventional Legendre 
transformation of thermodynamics : 

%-\-'^ (10) 

where T is the absolute temperature. This elegance is lost by Tsallis- 1 version which gives : 

^ = 1 = Z^-^ki5 (11) 

dU T ^ V ; 

where U is no longer a simple function dependent only on Z just like in the BGS caseQ. 

The Curado-Tsallis formalism once again gives Eq.(|^) with a change (g — 1) -^ (1 — g), and 
is successfully applied to systems showing non gaussian distribution laws as soon as its proposal. 
Nevertheless, the anomalous relation between the normalization Eq.(H) and the averages Eq.(P]) 



^In Tsallis — 1 version, we have U = — ^Wg X/ [^ ~ (? ~ l)/5ei]|' ^ , which nevertheless tends to BGS relation U — —-MglnZ 
when g — > 1. 



remains an open question. Rigorously speaking, supposing Eq.(|D, we logically write Eq.(|D for an 
observable average. In other words, Eq.(|^) means that pf is the observable probability but is not 
normalized and that Pi is an imaginary probability which is normalized but never used in practice, 
i.e. Pi is not observable if Eq.(^) holds. 
This second formalism, according to some authors P,|5|,pl|, shows some fundamental difficulties. 



For examples : 1) the average of a constant is not constant; 2) the total energy of two systems 
without any interaction is not the sum of the energies of each system; 3) the zeroth law does not 
hold; 4)the invariance of probability with uniform energy translation is missing. 

The third or last version of TDF is that mentioned at the beginning of this paper. It replaces the 
unnormalized average by a normalized one Eq.(RI) and so resolved the above mentioned problems. 



When writing this manuscript, we saw another proposal by Martinez et al |]13l concerning the 
energy constraint for entropy maximization. They propose using 

w 

Y.PKe^-U), (12) 

i 

as the constraint but with U given by Eq. (^ . This approach overcomes a mathematical difficulty of 
the third version of TDF, that is, the auxiliary function A{p) for entropy maximization : 

A = - ,_f' - «(Ep^ - 1) - Pi^f^ - U) (13) 

y^ i 2^i Pi 

does not necessarily have a maximum. Because, first, -;f-j- is not diagonal, i.e., J^ ^ ^ when 
i i^ 3- Second, even if we calculate only ^, we get : 

^ = -qp.{Z^-' - g/3pr^Z2«-2(e, - f/)]. (14) 

To ensure 44 < for a maximum entropy, we have to put \Z^^'-^ — ql3pj^~ Z^''^^(ej — U)] > 0, which, 
evidently, can not be always valid. Martinez et al overcame this difficulty with Eq.(|l2D leading to 



for i ^ j and -rr = —QPiZ^ '^ < 0. The entropy maximum is ensured. 



dpidpj I J d-pl 

This maximization still leads to Eq.(|l]) as probability distribution. In fact, this proposal is a 
composition of the Curado-Tsallis maximization with Y^ vi^i ^^nd another constraint Y^ vV^ ■ The 
latter can be called invariance constraint since its only role is to ensure the invariance of the resulted 
distribution with respect to uniform translation of the energy levels ti. Due to the same distribution 
functions, this scenario has the same essential characteristics as the last version of TDF, about 
which we will discuss some open questions in the following section. 

III. SOME QUESTIONS ABOUT TDF 

1. Let us begin by presenting an observation about the maximum entropy principle of Jaynes. 
This principle asserts that, in order to obtain the correct probability distribution, it suffices to 
introduce into the entropy maximization physical conditions or constraints related to observable 
quantities, that is, the normalization of the probability, the expectation, variance or higher 
moments. Supposed Eq.(^) as normalization, the only observable expectation value or higher 
moments must be defined with Eq.(^. The introduction of Eq.(|^) or Eq.(|^) does not conforms 
with this principle, because we do not know whether these averages really represent observable 



quantities or not. One may say that the successful apphcations of these formahsms confirm the 
observabihty of these averages. But in this case, Pi becomes non observable in turn. So the 
normalization Eq.(^ should disappear in the entropy maximization. This observability question 
deserves to be carefully studied. Indeed, if necessary, we could reject any old principle and 
introduce new ones under the condition that the resulted theory be useful and self-consistent. 

2. The second question was discussed by Raggio [|I^. If the expectation Eq. (^ satisfies the 
constraint associated with the linearity in the observables, i.e. 

X + y = x + y (15) 

for two independent observables x and y, it does not satisfy the linearity in the state (or 
distribution), i.e. 

x[Xp{l) + (1 - X)p{2)] = Xx[p{l)] + (1 - X)x[p{2)] (16) 

where p{l) and p{2) are two normalized probability distributions, nor does Eq. (P) which violates 
even Eq.(^). The origin of these violations is logically the discrepancy between the average 
Eq- (0j (or ^) and the normalization Eq. (^ according to which the joint probability 
is p{x + y) = p{x)p{y) and the probability summation is p{l + 2) = p{l) + p{2). If the 
average is calculated with p"^ instead of p, the problem is evident since p'^ is not normalized and 
p5(l + 2) 7^ pi{l) + pi{2). Indeed, 

-rx (^^^n x^ io^^ EaAp,(l) + (1 - A)p,(2)]%, 

xiXpiD + (1 - A)p(2)] = ^^[,^^(i) + (i_,)^^(2)]. ^''^ 



Eqs. ([T5| ) and (|16D are to be satisfied for that x, y, p{\) and p(2) be observable and have physical 
signification. This problem naturally leads us to the following question : does the observability 
of the normalization Eq. (§) is incompatible with that of the expectation value given by Eq. 
(I) or (i? 

3. The third question we would like to discuss is about the absence of analytic energy correlation. 
We know that the correlations in entropy can be calculated a priori with mathematical rigor. 
But the correlation in any other observable quantities can not be calculated with the average 
Eq. (I). 

Let us suppose an isolated system C composed of two subsystems A and B of which the 
distributions satisfy 

p,^{C)=p,{A)p,{B). (18) 

In nonextensive statistics, this hypothesis of multiplication law means that A and B are cor- 
related and gives the correlation term associated with energy. From Eq. (|T|), (^ and (0), we 
straightforwardly obtain : 

e.,(A + 5) - U{A + B) = [e.{A) - U{A)] Y.pI{B) + [e,{B) - U{B)] Y^p]{A) (19) 

+ {q-l)P[eU)-U{A)][e,{B)-U{B)]. 



Without additional hypothesis, this equahty does not lead to any explicit relation between 
the hamiltonians H{A + B) and H{A) or H{B) neither for the micro-state values Cj nor for 
the average U. Tsallis and coworkers 0,0 proposed neglecting the correlation term of any 
observable between the subsystems. So one can write : 

eij{A + B) = e,iA) + ej{B) (20) 

and 

U{A + B) = U{A) + U{B). (21) 

On the basis of this extensive energy approximation, the zeroth law of thermodynamics is 
claimed to be established for the foundation of nonextensive thermodynamics [|15],|16|. This 
problem will be discussed below 

Our question about Eqs. ( pOD and (^) is : if in general the correlation (nonextensive) terms of 
whatever observable or interactions can be neglected, what is the origin of the nonextensivity 
of entropy? We know that entropy should be a continuous function of the distributions which 
in turn are continuous functions of the observables. 

In addition, allowing Eqs. (pOD and (pT]), we lose the equalities Eqs.(|^) and Eg . ([l^) , which are 
crucial for the nonextensive theory. If Eq.(|l8]) fails, we cannot in fact find even the entropy 
correlation given by 

SiA + B) = SiA) + SiB) + ^^SiA)SiB). (22) 

k 

4. Zeroth law of thermodynamics can not be established without neglecting the correlation en- 
ergy. That is to say that the Lagrange multiplier (3 is no longer a "meter" measuring the 
thermodynamic equilibrium. We will show why. 

We take once again the above mentioned isolated system C composed of two subsystems A and 
B in equilibrium. From Eq. (^), we get, for a small variation of the total entropy : 

SS{A + B) = [l + ^i^SiB)]5S{A) + [1 + 1^5(A)]55(5) (23) 

Because 6S{A + B) = 0, we get 

Now, in order to proceed, we need the relation between U{A + B) and U{A) or U{B) but 
it does not exist. As mentioned above, U{A + B) can not be expressed in U{A) and U{B). 
So no relation can be found between 6U{A) and 6U{B) and, as a consequence, the derivative 
II can not be calculated from 6S{A + B) = 0. So supposing U{A + B) = U{A) + U{B) or 
5tJ{A) = —6U{B) proposed by Abe and Martinez et al |1T5| , [16| , we obtain 



or 

Z'''\A)f3{A) = Z''-\B)(3{B) (26) 

which is the generahzed zeroth law with Z'^~^(3 instead oi f3 = 1/kT as the measure of the 
equihbrium. Due to this approximate zeroth law, the implicit distribution function Eq. (P 
becomes explicit for systems in thermal equilibrium since it can be recast as 

[1 - (1 - g)A(e.. - U)]f^ 
Pi = ^ — (27) 

where A = f3Z'^~^ is now an independent thermodynamic variable. So the theory is reconciled 
with the old notion of thermal equilibrium to the detriment of the correlation between the 
components of the system. In our opinion, two important things are lost in this treatment : a) 
the generality of the nonextensive theory which should be able to tackle correlated problems 
in taking into account the interactions; b) the generality of the zeroth law which should hold 
within a theory without any condition, or all thermodynamic laws will become approximate. 



As a matter of fact, as mentioned above and argued by Guerberoff et al [|Tl|, with additive 
energy or non interacting subsystems, Eq.([T8|) and Eq.(|2^) do not hold so that the zeroth law 
Eq. ( p6D can not be established. Very recently, a new point of view about thermal equilibrium 
and nonextensivity shows that Eq.([T8|) or Eq.(^) is required by the existence of thermal equi- 
librium with Tsallis entropy and should be regarded as a basic assumption of the statistics for 
equilibrium nonextensive systems. So that appropriate energy nonextensivity satisfying Eq. (|18]) 
is absolutely necessary for the validity of zeroth law within Tsallis statistics [^,0. We will 
come back to this issue later in this paper. 

5. Now we discuss a mathematical problem. From Eq. (|1]), we can write : 

W 1 

i 

= E[l-(l-?)Me.-t/)]p' 

i 

w 

= z''$:pni-(i-9)A(e.-t/)]t 

i 

Considering the average defined in Eq. (0), Eq. (p8|) can be recast as : 

W q 

Z = Y.\^-{l-q)\{e,-U)Yf^ (29) 



or 



E[l - (1 - q)Ke^ - t/)]r = E[l - (1 - 9)Me. - U)]]-\ (30) 

i i 

This equality is a basic relation of the theory and must hold for arbitrary value of g, w, (3 and 
Cj. Now let us apply it to calculate the inverse temperature /3. 

From Eq. (pQ]), we can write 



Ep' = Z'-' (31) 

i 

and 

S = k 

1-q 



Z^-'i(U) - 1 
S = k^^ , (32) 



Then we calculate the following derivative : 

dS k dZ 
dU~^dU' 

First we take the Z given by Eq. (^, we obtain : 



(33) 



1 Q 

^ = kXZ^"^ = kp. (34) 



dU 
But if we take the Z of Eq. (^9]) , we obtain 

§ = '-^t[^-{l-^)He.-U)]r (35) 

which means 

W 29-1 

Z = g^[l-(l-g)A(e,-f/)]^^-V (36) 

i 

If we put Eq. (|3^) into Eq. (Q) and continue in this way for n times, we will find 

W q 

Z = Y.\^-{l-q)X{e,-U)]}-^ (37) 

i 

W 2q_ 

= ?E[l-(l-?)Me.-f/)]t" 

i 

w 

= q(2q-l)J2[l-il-q)\ie,~U)]:^- 

i 

w 

= q{2q - l)(3g - 2) ^[1 - (1 - g)A(e, - U)]^' 



2q~l 



3q-2 
-1 



4q-3 



{n+l)q — n 



= q{2q - l)(3g - 2)...(ng - n + 1) ^[1 - (1 - g)A(e, - U)]^ '~ 

i 

with n = 0, 1,2.... We create in this way a series of equalities which seem not to hold. For 
example, if we take the second equality of Eq. (^) and let g — > 0, the right-hand side will tend 
to zero and the left-hand side to Y^^ 1 = w. The result is lu ^ 0. This same result can also 
be obtained for g — >■ 1 if we take the third equality of Eq. (|37D and for g — > | with the forth 
equality and so on. These singular points in q value do not conform with the hypothesis that 
Eq. (^O]) is a general relation of the theory. It seems to us that these equalities are valid only 
when g — ^ 1 and Z becomes the BGS partition function. 



We can also study the equality Eq.(pOD in another way. Let us suppose that x is a continuous 
variable within < x < oo. So Z sometimes can be given by 

Z= [l-(l-g)A(a;-t/)]prfx. (38) 



or 







[l-{l-q)X{x-U)]^-''dx. (39) 



In this case, we should put g > 1 for Z to be calculated when x is large. The integration of 
Eq . (pSD is always finite. But Eq . (pOf ) needs g < 2 to be finite. If g > 2, the Z of Eq.(p8D can be 
calculated while that of Eq.(|39D diverges. This paradox naturally disappears for q —>■ 1. 



6. Another problem concerning the relation between the average U and the micro-state value e, 
or Z arises due to Eq. (pO]) . Usually, the U — Z relation (U = —-^InZ in BGS) can be found 
by introducing the distribution function Eq.(|I|) into the average calculus Eq.(^: 

U = ^^ (40) 

i 

= -|{|^ E[l - (1 - ^)Me. - U)]p - E f^[l - (1 - ?)A(e. - f/)]r } 

1 dZ „ 

which leads to, instead of the expected U — Z relation, 

as well as 

dS d Z^-'i-l \ dZ ^ 

ax = aA*^^r = '^^ = °- ''" 

So the micro- macro relation is impossible to be found if no mechanical quantity can be calculated 
from its microstate values. 

In addition, Eq.(^) gives rise to another problem similar to that discussed in the precedent 



part. Eq.(41) can be easily verified if we take the standard Z given by Eq.(|^). But If we take 



the Z in Eq.(|20), we get the following relation 



2q~l 



E(e. - f/)[l - (1 - g)A(e, - f/)]/- = (43) 



or 






If we repeat the same reasoning with the Z of Eg. (|36|) , we get 



U = "^^^^Pt- (45) 

2-^i Pi 



We can continue in this way with Eq.(40) and obtain : 



^ = ^^ (46) 

Z^i Pi 

EW 2(7—1 

^ i (^iPi 

Y^w 2q-l 
Z^i fi 



Ere^p 



3q-2 



y^w „3(3-2 
Z-ii fi 

ET e.pt'-' 



4g-3 

Ere.p: 

Ew nq-n+1 
i Pi 



nq—n+l 



which means ETi^i — U) = or U = Ef ^i/w if g = ^^ with n = 1,2,3.... i.e., we are led 
to the microcanonical case. On the other hand, if we calculate U from its definition Eq.(|]) for 
q = 1/2, we will be led to 

i i i 

or 

U=iT.e^-^T.e^)/iZ-^^), (48) 

i i 

which does not seem to be a microcanonical case. 

These mathematical problems discussed above seem to be related directly to the fact that the 
TDF Eq.(0) depends only on the difference Cj — U. So Eq.(0) can be regarded as the relative 
probability distribution with respect to the average energy of the system, but not the real 
probability with respect to a suitable zero-energy we must choose due to the energy translation 
variance of TDF. We will come back to this problem at the end of the present paper. We also 



wonder if the g-independent results for generalized ideal gas [13| have something to do with 
this relative probability combined with extensive energy approximation for "non interacting 
particles" . 

We have discussed some questionable points in the last version of TDF. Some other discussions 
concerning the relation between TDF of the third version and the law of large number are given in 
reference [|1^. Of course, these are questions that we have to study carefully. For the moment, we do 

10 



not see how to get out of these theoretical difficulties if we stay in the formalism with the conventional 
normalization and the expectation value Eq.(^. In what follows, we will show a possible way out. 
The main idea is to introduce into Tsallis theory the hypothesis of incompleteness of our knowledge 
(information) about nonextensive systems. 

IV. SOME CONSIDERATIONS CONCERNING STATISTICS 

Very recently, a possible alternative theoretical basis for TDF was proposed PJ20|]. The new 



formalism is based on a reflection about the conditions of physical application of the standard 
probability theory which is sometimes referred to as Kolmogorov probability theory ||2l[|. This 
reflection leads to a mathematically simpler and coherent nonextensive framework being capable of 
avoiding the problems discussed above. 

This new theoretical framework is referred to as incomplete statistics (IS). The basic idea of IS 
is that our knowledge about the states (their position and total number in the phase space) of a 
system may be incomplete and non exact. This is true at least for complex systems with unknown 
space-time correlations which can not be exactly described with analytic methods. So the equation 
of motion must be incomplete in the sense that some interactions are missing and its solution can 
not yield complete knowledge about the systems. On the other hand, Kolmogorov probability theory 
is founded on the hypothesis that we know all the possible states of the studied system or that the 
maximal information is complete and can provide definite answers for all questions that can be asked 
of the system. When we carry out a summation or an integration of probability in phase-space, we 
suppose that this is done over the possible states which can be determined by the equation of motion. 
This assumption is logical if and only if all interactions and space-time correlations are well-known 
or their unknown parts are negligible, as in the case of BGS or of other successful probabilistic 
sciences. 

V. INCOMPLETE NORMALIZATION 

If the incompleteness of our knowledge is not negligible, it is no longer sure that we sum over the 
possible states of the system simply because we do not know all of them. What we can do is to take 
the known states or events suggested by the equation of motion or by our knowledge. Their number, 
say f , may be greater or smaller than w, the real number of all possible states. As a consequence, 
the normalization condition is reduced to 

jlp^ = Q^l (49) 

where Pi is the true probability which can not sum to one. The necessity to introduce this nonex- 
tensive or nonadditive probability is first noted by economists |jl[]. This probability in Eq.(^9[) was 
referred to by Renyi as incomplete probability distribution |^ because the values of the random 
variable of this distribution do not constitute a complete (exhaustive) ensemble (i.e. v < w). 

It should be noticed that Q is a constant depending on the studied system. Now in order to apply 
conventional probability theory, Eq.(^9|) has to be renormalized to get a calculable "probability" 
related to pi as well as to Q representing the nature of the system. 

Our proposal is to write 

t% = tpl (50) 

11 



which means 

Ep? = i,(ge[0,oo]) (51) 

and, for internal energy, 

U = J2ptei. (52) 



j=i 



Since Pi is the true probabihty and satisfies < pj < 1, we have to set < g < cxd. q = should 
be avoided because it leads to pj = for all states. We see from Eg . (^Of) that Q = 1, Q < 1 and 
Q > 1 means g = 1, g < 1 and g > 1, respectively, g is directly related to Q (e.g., for microcanonical 



■u— In 

help to understand empirical g values for nonextensive systems. 



distribution, g = ^^ j^"^ „ ) and, in this way, to the unknown correlations or information. This may 



VI. NONEXTENSIVE STATISTICS 

In this section, we present a method based on the incomplete normalization Eq.(|5T|) in order to 
get TDF. 

A. Nonextensive information and entropy 

On the basis of the hypothesis of the nonextensivity (nonadditivity) of entropy or other quanti- 
ties, we proposed using a generalized logarithm function as a generalized Hartley formula for the 
information measure I(N) required to specify one element in a system containing N elements : 

N^ -1 

I{N) = (53) 

9 

where 5^ is a real number. This means the following nonadditivity : 

/(A^i X N2) = /(A^i) + I{N2) + gI{N^) x /(A^a). (54) 

where I{Ni x A^2) is the information needed to specify simultaneously 2 elements, each being in a 
subsystem 1 or 2. We see that the parameter g is a measure of nonextensivity. If (? — ^ 0, I{N) -^ InN 
and the information becomes extensive (additive). Eq.(|5^), Eq.(p3D or Eg . ([5^) , plus the other axioms 



used by Shannon pj2^, lead to a nonextensive entropy 

S = -kY^pp = -k^^i^l^ -. (55) 

f^i 9 9 

We should ask that the above entropy recover the Gibbs-Shannon one 5" = —kY^^=iPilnpi for 
g = 1. This constraint on nonextensive entropy is logical because g = 1 or Q = 1 implies a complete 
knowledge about the studied system and short range interactions or correlations. In this case, we do 
not have any reason for holding the nonextensivity. So (7 = when g = 1 and g should be monotonic 
function of g, which ensures the monotonic g-dependence of entropy. 
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B. generalized distribution function 

For maximum entropy, we write the following auxiliary function : 



Let ^ = 0, we obtain 

dpi ' 



A = ^^ + «7Ep? - «/5'EP?e.. (56) 



P^ = ^^^ (57) 



with 



^' = E[l--e.]r. (58) 

i ' 

Now it should be asked that the distribution Eq.(^) become the BGS exponential distribution for 
g = 1 or g{q) = 0, which means that Eq.([57|) should be a generalized exponential function, that is 

?MM^ . _^'e.. (59) 

9 

This straightforwardly leads to 7 = 1/ g{q) and 

Pi = -^ — (60) 

From Eq.(^), it can be shown that 

^ = -[9iq) + q]pf'^-''-'- (61) 

If we want that the distribution Eq.(^) be a maximum entropy (or minimum energy ^^) distribu- 
tion, we should put 

g{q) + q>0 (62) 

which ensures y^ < for a maximum entropy with Eq.(^). This means that the curve of g{q) is 



dpf 

situated above the straight line of g{q) = — g. If we impose the monotonicity of S and g{q) with 
respect to q, we also have 

^ < 0. (63) 

dq 

Eq.(|60D is the generalized nonextensive distribution function. Comparing it to Eq.(|I]), we see that 
TDF corresponds to g{q) = 1 — q. This choice is the simplest one that satisfies the two conditions 
Eqs.(p^) and ( pBD and also yields the simplest nonextensive entropy : 

S = -kY^pf-^ = -A:E^^^ = k ^^'='^' . (64) 

U 1-9 f^i 1-9 1-9 

13 



with the simplest generahzed distribution function 

which is the TDF in IS version. It should be remembered that a different forms of g{q) may lead 
to different distributions. 

Note that the distribution Eq.(^) is in the same form as the TDF in Curado-Tsallis version but, 
due to the incomplete normalization, has a different partition function given by 

Z« = ^[l-(l-g)/3'e,]^^. (66) 

i 

The consequences of this change will be discussed in the following sections. 



We would like to mentioned here that Kaniadakis p4| , |25| proposed a new generalization of BGS 
(^-deformed statistics) on the basis of nonlinear kinetics in low density gas systems with normalized 
distribution function /. It is argued that the g-variance of the g-exponential function of TDF 
suggests to write f = pi and so Eqs. (|5T|) and (|52D . This work gives Eq.(^) as single particle 
distribution from the nonlinear kinetics theory. This result is in accordance with the conclusion of 



reference |T^ that TDF is an exact distribution for both many-body system (nonextensive) and 



correlated single body according to Eq.(|T8D prescribed by thermodynamic equilibrium |[T7|JT8[| . In 
this sense, the hypothesis of low density gas is no longer necessary, and the K-deformed statistics can 
be regarded as a valid theory for any nonextensive gas system in equilibrium having Tsallis entropy. 

C. Generalized distribution and nonextensivity 

We consider again the total system C composed of two subsystems A et B in interaction. By 
Pij{C) we denote the probability that C is at the product state ij while A is at the state i with a 
probability Pi{A) and B at j with a conditional probability Pij{B \ A). We can write 

p,,{C)=p,iA)p,^iB\A), (67) 

or 

pUC)=pKA)pUB\A), (68) 

which means 

[l-giq)(3'e,,{C)]f^ ^ [1 - g{q)P'e.{A)]f^ [1 - g{q)(3'e,,mF (.n) 

Z'i{C) Z<i{A) Zi{B I A) ^ ' 

e,,{C) = e,{A) + e,,{B) - g{q)(3'e,{A)e,,{B). (70) 

UiC) = UiA) + U{B) - g{q)(3'UiA)U{B). (71) 



or 



and, from Eq. (|5^) , 



Since S is an observable as the others, Eq . ([71]) must hold for S as well. Indeed, if we put the 
multiplication law Eq.(|67D into the entropy Eq.(^), we get 
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S{C) = -k ^'^P'^ \^^ (72) 



9{(l) 
-k 



with 



9{q) 9{q) 

k 9{q) 9{q) 

SiA) + S,{B I A) - ^S{A)S,{B \ A). 

k 



S,B I ^) ^ _,S.^r'"(^M)-\ p3) 



Following the idea of Abe et al |2^, we can define a conditional entropy of B as follows 

S{B\A) = Y.p^{A)S,,{B\A). (74) 



From Eq.(^), we simply obtain 



1 - SMS(A) 
or 



S(BlA) = ZPnA)f^^^^ (75) 

_ S{C) - S{A) 



S(C) = S(A) + S(B I A) - ^S{A)S{B \ A), (76) 

k 



which has exactly the same form as Eg. ([7TD . It should be noticed that this equation means that 
S{B I A) = Si{B I A), i.e. the conditional entropy of B associated with a microstate i of A is in 
fact independent of i. It is a little surprising to see this property with two correlated subsystems. 
But if we look back to the origin, we find that this is simply the consequence of the postulate 
Eq. ( p^ which supposes a nonextensive term in the form of product of two sub- informations. This 
is naturally a special choice of the nonextensivity described by the generalized Hartley formula. 
If we choose a different nonextensive information measure, the things will be different. Recently, 
the proposition of an entropy pseudoadditivity |T^ prescribed by thermal equilibrium shed light on 
this problem. We understand that, supposed the generalized Hartley formula (or Tsallis entropy), 
Eq.(^4|) and so Eq. (|2^) are prescribed by the existence of thermal equilibrium. In other words, 
without these relations, a composite system with correlated (not independent!) subsystems can not 
have stable equilibrium state. So any exact discussion about equilibrium systems must conform with 
the factorization of compound probability Eq . (|I8|) . As a consequence, the idea of additive energy 
with neglected correlations [l^Jl^ or non interacting systems [jll|, incompatible with the spirit of 
nonextensive statistics, becomes unnecessary and should be rejected. The definition of temperature 
has to be revisited on the basis of nonextensive energy satisfying Eq . (p!8D . It is what we are doing 
in the following section. 
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VII. THERMODYNAMIC RELATIONS 

To give statistical interpretation of thermodynamics, a well defined temperature related to stable 
thermodynamic equilibrium (maximum entropy or minimum energy) is needed. In this section, we 
will present briefly some consequences of IS. All the discussions are based on Eq. ([TsD and compatible 
entropy and energy pseudoadditivities. 

A. Zeroth law and generalized temperature 

First, what is /?' in the distribution Eq.(|60|)? In BGS, /?' = /?= ^ is the inverse temperature and 
the first law of thermodynamics can be written as 

dU = dQ + W (77) 

or 

dU = TdS + YdX (78) 

where W is the work done by Y, a generalized exterior force (e.g. pressure -P), X the correspondent 
displacement (e.g. volume V), and dS = ^ the thermodynamic definition of entropy. When X 
remains constant, we have dU = dQ and 

On the other hand, within IS, from Eqs.(|55|) and (|60D, we obtain 

Z-9M -1 , . 

S = -k — + kp'Z-a^'^^U. (80) 

9{q) 



or 



1 - Z^^'^ 

S' = -k + kp'U. (81) 

l-q 

with S' = Z9(i)S = Z^-'iS. This leads to 

W-'^' (82) 

Now we are showing that (3' still measures thermal equilibrium at maximal entropy or minimal 
energy. Let us take the nonextensivity relation Eq.(^) and calculate a small variance of energy 
U{C): 

5U{C) = [1 - g{q)P'U{B)]5U{A) + [1 - g{q)^'U{A)]5U{B). (83) 

At equilibrium or energy minimum, 6U{C) = 0, we obtain 

[1 - g{q)P'U{B)]6U{A) = -[1- g{q)f3'U{A)]SU{B). (84) 

Putting this equation into the entropy maximum relation Eq.(|^) [in which {1 — q) should be replaced 



by {q — 1) due to the IS version of entropy nonextensivity Eq.(|76|)], we get 
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1 - g{q)P'U{A) dSjA) ^ l-g{q)(3V{B)dS{B) 
l-9MS{A) dU{A) i-SMs{B) dU{By ^ > 

With the help of Eg. (pOD , we can estabhsh 



or 

dS'{A) dS'{B) 



87) 



dU{A) dU{B) 

which means 

(3'{A)=(3'{B). (88) 

This result can also be obtained in another way. Multiplying Eq.([76D by Z^^'^[C) and considering 
Z{C) = Z{A)Z{B), we obtain 

S'{C) = Z^-\B)S'{A) + Z^-''{A)S'{B) - ^S'{A)S'{B). (89) 

On the other hand, from Eq.(^), it is straightforwardly to verify that, for maximum entropy at 
equilibrium, 6Z{C) = and 6S'{C) = 0. This leads Eq.(|P|) to 

[Z^^'iB) - ^^S'{B)]6S'{A) = -[Z^-'{A) - ^^S'{A)]6S'{B). (90) 

Now with the help of Eqs.(|3) and (|8|), we obtain Eq.(|3) and Eq.(|8D . 
So /?' remains the meter of stable equilibrium state. We can define a generalized temperature 

where T' = T if g = 1 or g{q) = 0. 

B. Some other relations 

Considering Eq.(p2D and the energy conservation law Eq. (|77D , we see that dS' should be the 
measure of heat transfer. We have to write 

dQ = ^, (92) 

or 

dQ = T'dS'. (93) 

Now the first law of thermodynamics Eq. ( [78D should be written as follows 

dU = T'dS' + YdX. (94) 

The free-energy F have to be defined as 
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dF = -S'dV + YdX (95) 

or 

F = U - T'S', (96) 

which leads to, with the help of Eq. (^) : 

F = -kT . (97) 

Considering Z{C) = Z{A)Z{B), we easily obtain : 

FiC) = F{A) + F{B) - ^-^F{A)F{B). (98) 



We also have, for the heat capacity 



^ _dQ_ dS' _ d'F 



and for the generalized force 

Y = i 

dX 



dF 
y = {^}t'. (100) 



VIII. THE FUNDAMENTAL PROBLEMS REVISITED WITHIN INCOMPLETE STATISTICS 



What about the problems discussed in section |T| if we consider I SI Let us examine them one by 
one. 

1. The first problem of the incompatibility between normalization and expectation value does not 
exist any more because the incomplete normalization Eq. (|51|) is compatible with the expectation 
value Eq.(|2D. 



2. The second problem was that the expectation value Eq.(^ was not linear in the distributions 



shown by Eq.(0). This problem can be avoided within IS thanks to Eqs.(0) and (52) : 



x[Ap(l) + (1 - A)p(2)] ^ ^^[,,.(i) + (i_,),.(2)] (10^) 

i i 

= Ax[p(l)] + (1 - X)x\p{2)] 

3. The problem of the absence of analytic correlation in energy and other observable quantities 
can be resolved by Eqs. (l70| ) and ([7l|) , which guarantee a nonextensive statistical theory with 
mathematical rigor. 

4. The problem of the zeroth law of thermodynamics was resolved in the above section. 

5. The fifth problem discussed in section |T| is due to the peculiar equality Eq. ( ^OD which in turn 
is due to the average defined in Eq. (^. In IS, this problem does not exist any more thanks 
to the incomplete normalization and the concomitant expectation value Eqs.(|5TD and (0). 
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6. In IS, the U — Z relation can be found from Eqs. (|52D and ([65| ) : 



U = J2p'e^ (102) 



1 ^ 

1 9 



^e,[l-(l-g)/5'e,]^^- 



1 dZ 

where 

Z = ^[l-(l-g)/3'e.][^. (103) 

i 
IX. ENERGY INVARIANCE OF THE DISTRIBUTION 



A crucial problem of the nonextensive distribution Eq. ( pq) is that it is not invariant with respect 
to uniform translation of energy spectra Cj. If we replace Cj by e, + C where C is constant, Eq.(|65D 
becomes : 

Pi = ^ — (104) 

with 

Z^ = ^[1 - (1 - g)/?'(e. + C)]f~^ (105) 



which is not same as Eq. (|65D , excepted that q = 1. This problem worries enormously scientists 
Q. It is known that thermostatistics takes into account only energies relative to thermodynamic 
movements and that the choice of energy-zero is never a problem in BGS because the theory has 
exponential distribution and is invariant with uniform energy translation. But with TDF, Eq.( |104[ ) 
implies that the properties of a gas may depend on the translation speed or on the location of the 
container. Although in practice we can always choose the usual energy-zero as for BGS, this peculiar 
theoretical property of TDF is somewhat unusual and disturbing. Avoiding this puzzling variance 
of distribution has been one of the motivations of the third version of Tsallis theory. Now the 
following questions arise : Why do we have to fix only one choice of zero-energy in TDF to avoid 
the container-dependence of the statistics? Why does the nonextensivity lead to this theoretical 
property? Is it really something to be avoided? 

As is well known, BGS is an extensive theory which holds only for systems with weak or short 
range interactions. In addition, the invariance of BGS is based on the classical mechanics which 
leaves the interaction potential completely arbitrary. But TDF is a theory for solving problems with 
complex interaction or correlations we often ignore. So we are not obliged or it is not advisable to 
identify the variance of TDF to that of BGS. The energy translation invariance is not an universal 
characteristic of physical theories. As a matter of fact, this property of arbitrary potential energy 
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disappears even in classical mechanics and BGS if we consider the relativistic effect. We can not 
add constant into the energy E oi a, system because E is related to its total mass M by 



E = Mc^ 



(106) 



where c is the light speed. E is not arbitrary because M can not be changed arbitrarily. 

If we accept the the variance of TDF, we have to choose a definite zero potential energy. Let us 
suppose a non negligible interaction or correlation between all the elements of mass rrii and energy 
Cj (z = 1, 2, 3, ...) of a system. Let V be the total interaction energy. M can be given by : 



M = Y,m + V/c^ 



(107) 



or 



Mc^ = Y^ruic^ + V = Y.ei + V 



(108) 



with Cj = rriiC for i element. It is obvious that neither M, nor Cj and V can be changed. If not, 



the variance of the theory would be perturbed |2^. According to this discussion, a possible choice 
of zero potential energy corresponds to the following case 



M = E 



rrii 



(109) 



or ^ = 0. This condition may correspond in some cases to infinite distance between the elements 
of the system (e.g. an atom when we consider the internal energy between the electrons and the 
nucleus), and in other cases, to special positions of this elements (e.g. equilibrium position of the 
atoms in crystal lattice if we are interested in their vibration). This is in fact just what we usually 
choose with BGS. 

It is worth emphasizing that some of the problems discussed in section III seem to be related to the 
distribution functions invariant through energy translation. Because, if we apply the maximization 
method of Martinez et al |T^, i.e. to introduce J2T Pi^ ^s the invariance constraint into the auxiliary 
function Eq.(p3), we obtain : 



Pi 



[1 - (1 - g)/3'(e. - U)];- 



with 



Z'^ = Y.l^-il-q)P'{e,-U)]^- 



(110) 



'111^ 



which is invariant through energy translation and different from Eq.(|T|) only by the partition function. 
We easily find that the problems 3 to 6 take place again with questionable equalities similar to 
Eqs.([T9|), ( p4D and (^) which make it difficult to establish the zeroth law without any approximation 
and lead to other puzzling equalities like Eqs.(P7|), (^D|) and (^Bj). Does the distribution invariance 
inevitably lead to the theoretical peculiarities? This seems an interesting topic which is beyond the 
range of the present work. 
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X. CONCLUSION 

We have shown some observations about the actual nonextensive statistical theory. The prob- 
lems discussed reveal that, with the conventional normalization Eq.(^ and the expectation value 
Eq.(^, the generalized exponential distribution, though very successful in many applications, cannot 
be obtained with convincing theoretical approach and so the nonextensive statistics shows peculiar 
properties which seem difficult to be avoided. We have shown that it was possible to overcome these 
difficulties if we introduced the concept of incomplete information with suitable normalization and 
expectation. This approach allows to establish TDF in a consistent way with only a different parti- 
tion function. New nonextensive thermodynamic relations were deduced on the basis of generalized 
definitions of heat and temperature. It is argued that the energy invariance should not be considered 
as a necessary property of TDF. The connection between the nonextensivity and the energy shift 
dependence of TDF remains to be understood. 
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